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Interpolation Newton’s Forward Difference Formula 
Let y =f (x) be a function which takes values fxg) fxgt h), fxgt2h), ..., corresponding 


to various equi-spaced values of x with spacing h, say xo Xo + h, Xot 2h, 
Suppose, we wish to evaluate the function f (x) for a value Xot ph, where p is any real 


number, then for any real number p, we have the operator E such that 
E” f(x) = f Œ+ ph). 
S + ph) = EP fa) = A+A) f) 


=| pat E-D a PP MPA) 99 | Fea) 
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F(X + PA) = f(x) + PAF (%) 

—] —1)\(p-2 
A e EOT p(p ae ) A fC) 
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n! 

This is known as Newton’s forward difference formula for interpolation, which gives the 
value of Ko + ph) in terms of Axo) and its leading differences. 


+: A" f (xo) + Error 


This formula is also known as Newton-Gregory forward difference interpolation formula. 
Here P=(X-Xp)/h. 


An alternate expression is 


(p-1) (p-1)(p-2) 
yantan A’y, + a A 
+ PRED ND Ary, + Error 
n.: 


If we retain (r + 1) terms, we obtain a polynomial of degree r agreeing with y, at 


x0, xl, ..., XT. 
This formula is mainly used for interpolating the values of y near the beginning of a set of 
tabular values and for extrapolating values of y, a short distance backward from Yo 


Example: 
Evaluate f (15), given the following table of values: 


x Ena Ree E eas eS 


¥=f@)| 46 


Solution: 
The forward differences are tabulated as 
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We have Newton’s forward difference interpolation formula 
P(p-l) Ka p(p-Ip-2) AN 


PERDA 31 


, Pie-Vp-nt) 
n! 


A” yo + Error 


Here we have 
x, =10, y, =46, Ay, = 20, 


Ay =-5, A’y,=2, A*y, =-3 


Let Yis be the value of y when x = 15, then 
BM _15-10 


p= =0.5 
h 10 
f(5) = y; = 46+ (0.5)(20) + Coe (-5) 
fi (0.5)(0.5 ee 2) (2) + (0.5)(0.5 “uw — 2)(0.5—3) (3) 


= 46+104+0.625+0.125+0.1172 
Ff (1S) = 56.8672 correct to four decimal places. 


Example 


Find Newton’s forward difference, interpolating polynomial for the following data: 


x (0.1 (0.2 (0.3 0.4 (0.5 


peccsccccescescccssecccoscocss: dp E EE PEE POE SOE 


y=f@) 14 156 1.76 2 


Solution; 
The forward difference table to the given data is 
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A A A EEE PAEAN EEEE EEEE EASA SAATA 


(oral o.1i4 
pei ae ess 
iP coke Ee i76 


rd th 
Since,3 and4 leading differences are zero, we have Newton’s forward difference 
interpolating formula as 


-1 
Y = yo + PAVo PO a2, 
In this problem, 
x, =0.1, ya =1.40, 
Ay, =0.16,A’y, = 0.04, 
and 


x—0.l lox] 


p = 
Substituting these values, 


y= f(x) =1.40+(10x-1)(0.16)+ 


TE E (0.04) 


This is the required Newton’s interpolating polynomial. 
Example 
Estimate the missing figure in the following table: 


Solution 


Since we are given four entries in the table, the function y = f (x) can be represented by a 
polynomial of degree three. 


A’ f (x) = Constant 
and A* f(x)=0, Vx 
In particular, 
A* f(x) =0 
Equivalently, 
(E-l)* f(a) =0 
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Expanding, we have 
(E* —4E° +6E° -4E +1) f(x,) =0 
That is, 
F(X) -4F 0G) F OF x) —4F (01) + FQ) = 0 


Using the values given in the table, we obtain 
32-4 f(x,)+6x7-4x5+2=0 


which gives f (x3), the missing value equal to 14. 


Example 

Consider the following table of values 

Xx :2, 3 A 5 6 
F(x) .2304 .2788 3222 .3617 .3979 
Find f (.36) using Newton’s Forward Difference Formula. 

Solution 


04| 03222] oo3ə5| -ooo33| | | 
05| osez] oos] J J | 
| 


| 06] osso] | Ù | 


4! n! 
Where 
% =0.2, y= 0.2304; Ay, =0.0484, rom _036-02_,, 
A?y, =—0.005, A®y,=0.0011,A‘y,=-0005" h& O1. 


y= 02304+1.60.0484) +S (0.005) eee Oo p SS DEL AA 0005 


= 0.2304 +.077441 -024+ LE 001 4 ALOE I™D (0995 


= 0.3078 —.0024 —.00007 —.00001 
=.3053 
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Example 

Find a cubic polynomial in x which takes on the values 

-3, 3, 11, 27, 57 and 107, when x = 0, 1, 2, 3, 4 and 5 respectively. 
Solution 

Here, the observations are given at equal intervals of unit width. 

To determine the required polynomial, we first construct the difference table 


Difference Table 
x Ta Cx] ar Cx] ; ve F(x) 
eee O ks Ge a eee 


th l ; ; 
Since the 4 and higher order differences are zero, the required Newton’s interpolation 
formula 


F(% + ph) = f(%) + DAF (%) + PPDA FC) + Pen at =) 3 £(x,) 


Here 
Js X—=X _ x-0 ats 
h 1 
Af (x) =6 
A f (%)) =2 
A? f (x)) = 6 
Substituting these values into the formula, we have 


Prange = =o =A) 


(6) 


f(x) = 2° -2x? + 7x-3, 
The required cubic polynomial. 
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